Traffic flow of interacting self-driven particles: 
rails and trails, vehicles and vesicles 
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One common feature of a vehicle, an ant and a kinesin motor is that they all convert chemical 
energy, derived from fuel or food, into mechanical energy required for their forward movement; 
such objects have been modelled in recent years as self-driven "particles". Cytoskeletal filaments, 
e.g., microtubules, form a "rail" network for intra-cellular transport of vesicular cargo by molecular 
motors like, for example, kinesins. Similarly, ants move along trails while vehicles move along lanes. 
Therefore, the traffic of vehicles and organisms as well as that of molecular motors can be modelled as 
systems of interacting self-driven particles; these are of current interest in non-equilibrium statistical 
mechanics. In this paper we point out the common features of these model systems and emphasize 
the crucial differences in their physical properties. 
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I. INTRODUCTION 
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In the recent years non-equilibrium statistical mechan- 
ics has found unusual applications in research on traffic 
flow of various different types of objects. In this pa- 
per we consider mainly three different examples of such 
traffic, namely, (a) vehicular traffic 0, |E 0, (b) ant- 
traffic on ant-trails 0, IE IE 1 an d ( c ) intra-cellular traf- 
fic of molecular motors carrying vesicular cargo moving 
along cytoskeletal filaments 0- Most of these models 
are essentially generalizations or extensions of the Asym- 
metric Simple Exclusion Process (ASEP) 9], which is, 
to our knowledge, the simplest model of systems con- 
sisting of interacting driven particles; the general aim of 
these investigations is to understand the interplay of self- 
organized structures and transport in systems driven far 
from equilibrium 0, 0, 0] . 

The aim of this article is (a) to summarize the main 
results of recent works on all the three systems men- 
tioned above, elucidating the nature of various types 
of quenched randomness, (b) to present the challenging 
open problems, and (c) to indicate the possible trends 
of future developments in this frontier area of interdisci- 
plinary research. 

The common modeling strategy is to represent the 
motile elements (i.e., vehicles, ants, molecular motors) 
by "self-propelled" particles which convert chemical en- 
ergy (derived from fuel or food) into the mechanical en- 
ergy required for the forward movement. In such generic 
models, the mutual influences of the motile elements on 
the movements of each other are captured by appropri- 
ate inter-particle interactions. In the spirit of the lattice 
gas models, the track for the traffic movement (i.e., the 
highway lane or ant trail or cytoskeletal filaments) are 
represented as discrete lattice of "cells" each of which 
can accomodate at most one partcle at a time. The dy- 
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FIG. 1: TASEP with periodic boundary conditions. The hop- 
ping probability of the particles is q. 



namical laws governing the forward movement of the self- 
propelled particles in such "particle-hopping" models arc 
usually formulated as "update rules" in terms of cellular 
automata (CA) 



II. VEHICULAR TRAFFIC 

The simplest model of interacting self-driven particles 
is the so-called totally asymmetric simple exclusion pro- 
cess (TASEP) HIIEII3- Imposition of open boundary 
conditions leads to richer physics as compared to those 
for the corresponding model with periodic boundary con- 
ditions. The states of the system are updated either in 
parallel or in a random-sequential manner following rules 
which will be explained later in this section. 



A. TASEP with periodic boundary conditions 

TASEP with periodic boundary conditions is sketched 
schematically in fig^ I n the original formulation of 
TASEP, the states of the system were updated in a ran- 
dom sequential manner where a particle is picked up ran- 
domly and moved forward by one lattice spacing, with 
the hopping probability q, provided the site in front of 
the particle is empty. 
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FIG. 2: Fundamental diagrams of the TASEP with peri- 
odic boundary conditions and (a) random-sequential updat- 
ing (solid curve) (b) parallel updating (dashed curve), both 
for q = 0.75. 



given by 



and 



F ^2 



F r = q c (1 — c) 



1 - \J\ - M c (1 - c) 



(1) 



(2) 



respectively; the two expressions Q and (J2J) , for identical 
hopping probability q — 0.75, are compared in figGl Note 
that both the expressions and J5J exhibit particle- 
hole symmetry, i.e., these are symmetric under the in- 
terchange of c and 1 — c. In section IIIII we shall show 
how this symmetry is broken in a model of ant-traffic on 
ant-trails. 
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FIG. 3: TASEP with open boundary conditions. 



The Nagel-Schreckenberg (NS) model [14j is a mini- 
mal CA model of vehicular traffic on idealized single-lane 
highways; the maximum possible (discrete) speed of the 
vehicles is V max . However, in the special case V rnax — 1 
this model reduces to the TASEP with parallel updating. 

The exact fundamental diagrams of the TASEP with 
random-sequential updating and parallel updating are 




FIG. 4: Phase diagram of the TASEP for parallel dynamics. 
The inserts show typical density profiles. 



B. TASEP with open boundary conditions 

The open boundary condition is, however, closer to 
the real vehicular traffic on a stretch of highway. If open 
boundary conditions are imposed on the TASEP, addi- 
tional rules must be specified to regulate the entry and 
exit of the particles at the two boundaries of the finite 
system. Usually, these are specified as follows: if the first 
site at the open point of entry is empty it is filled with 
probability a whereas particles occupying the last site at 
the point of exit hop out of the system with probability 
P (see figEJ). 

The open boundaries break the translational invari- 
ance of the system and give rise to stationary states 
with non-trivial density profiles. Such model systems 
have been investigated thoroughly over the last decade 
from the point of view of fundamental principles of non- 
equilibrium statistical mechanics. In contrast to equilib- 
rium systems with short-range interactions, such driven 
non-equilibrium systems can exhibit transitions from one 
dynamical phase to another, even in one-dimension with 
only short-range interactions, with the slight change of 
boundary conditions 0, HEj . 

The typical phase diagrams of the TASEP with open 
boundary conditions are sketched in fig01 the qualitative 
features of the phase diagram of TASEP is practically 
independent of the nature of the dynamics. In the low- 
density phase A the flux is independent of (3 and lim- 
ited only by a. On the other hand, in the high-density 
phase B the flux is independent of a and determined by 
(3. However, in the maximum flux phase C the current 
is independent of both a and f3. Moreover, both the 
high- and low-density phases can be subdivided into two 
phases each, namely, AI, All and BI, BII, respectively; 
these subphases are distinguished by the asymptotic be- 
haviour of the density profiles at the boundaries (see the 
insets in fig0J. In section Hvl we shall see how this phase 
diagram is modified in models of intra-cellular transport 
of vesicular cargo by molecular motors. 
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each stage the dynamical rules are applied in parallel to 
all ants and pheromones, respectively. 

Stage I: Motion of ants 

An ant in cell i that has an empty cell in front of it, i.e., 
Si(t) — 1 and Si-|-x(t) = 0, hops forward with 



probability 



Q if a l+1 (t) = 1, 
q if <J i+ i(t) = 0, 



(3) 



where, to be consistent with real ant-trails, we assume 
q<Q. 



FIG. 5: Schematic representation of typical configurations; 
it also illustrates the update procedure. Top: Configuration 
at time t, i.e. before stage I of the update. The non-vanishing 
hopping probabilities of the ants are also shown explicitly. 
Middle: Configuration after one possible realisation of stage 
I. Two ants have moved compared to the top part of the 
figure. Also indicated are the pheromones that may evaporate 
in stage // of the update scheme. Bottom: Configuration 
after one possible realization of stage II. Two pheromones 
have evaporated and one pheromone has been created due to 
the motion of an ant. 



III. ANT TRAFFIC ON TRAILS 

Ants communicate with each other by dropping a 
chemical (generically called pheromone) on the substrate 
as they crawl forward 0, 0|. Although we cannot 
smell it, the trail pheromone sticks to the substrate long 
enough for the other following sniffing ants to pick up its 
smell and follow the trail. In our recent papers 0,0 we 
have developed a particle-hopping model, formulated in 
terms of stochastic CA, which may be interpreted as a 
model of unidirectional ant-traffic on a trail. Rather than 
addressing the question of the emergence of the ant-trail, 
we have focussed attention on the traffic of ants on a trail 
which has already been formed. 

Each site of our one-dimensional ant-trail model rep- 
resents a cell that can accomodate at most one ant at a 
time (see Fig. |5J. The lattice sites are labelled by the 
index i (i — 1,2, ...,£); L being the length of the lat- 
tice. We associate two binary variables Si and o~i with 
each site i where Si takes the value or 1 depending on 
whether the cell is empty or occupied by an ant. Simi- 
larly, o~i = 1 if the cell i contains pheromone; otherwise, 
o~i = 0. Thus, in contrast to TASEP, we have two subsets 
of dynamical variables in this model, namely, {S'(t)} and 
{a(t)}. The instantaneous state (i.e., the configuration) 
of the system at any time is specified completely by the 
set ({$},{*}). 

Since a unidirectional motion is assumed, ants do not 
move backward. The forward-hopping probability of an 
ant is higher if it smells pheromone ahead of it. The 
state of the system is updated at each time step in two 
stages. In our ant-trail model with parallel dynamics, at 



Stage II: Evaporation of pheromones 

At each cell i occupied by an ant after stage I a 
pheromone will be created, i.e., 

o-<(t + l) = l if Si(t + 1) = 1. (4) 

On the other hand, any 'free' pheromone at a site i not 
occupied by an ant will evaporate with the probability / 
per unit time, i.e., if Si(t + 1) = 0, o~i(t) = 1, then 

. > _ J with probability /, , . 

cr^r -t- -l; < [1 with probability i _ f K°) 

We have also considered another version of our ant- 
trail model where the states of the system are updated 
in a random-sequential manner rather than in parallel 
Note that in both the cases, because of the periodic 
boundary conditions, the dynamics conserves the number 
N of ants, but not the number of pheromones. 

This model is related to several other models. For 
example, in the limits / — ► and — * oo this ant-trail 
model reduces to TASEP with the hopping probabilities 
q and Q, respectively. Moreover, the ant-trail model may 
be regarded as the opposite limit of the bus-route model 
UBEil (see ref.pl for the detailed comparison). Further- 
more, the ant-trail model also has some similarities with 
the par ticle- hopping models of human trails of pedestri- 
ans [uim. 

The typical fundamental diagrams of the ant-trail 
model |4j with parallel dynamics are shown in figlHI the 
corresponding results in the case of random-sequential 
updating are qualitatively similar ||. The unusual 
shapes of the curves observed over a range of / are con- 
sequences of the non-monotonic variation of the average 
speed of the ants with their density on the trail (see figCJ. 

Both the ordinary mean-field theory (MFT) , which ac- 
counts for the exact fundamental diagram of the TASEP 
with random-sequential updating, and 2-cluster MFT 
0,122, which succussfully predicted the exact fundamen- 
tal diagram of the TASEP with parallel updating, fail to 
capture even the qualitative features of the fundamental 
diagrams of the ant-trail model shown in figEl @ • How- 
ever, a heuristic MFT, described in ref.Q, captures at 
least the qualitative features of the observed flow prop- 
erties of our ant-trail model. 
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FIG. 6: Fundamental diagram of the ant-trail model with par- 
allel updating for the parameters Q — 0.75, q = 0.25. The dis- 
crete data points corresponding to / = 0.0005(0), 0.001(o), 
0.005W, 0.01(A), 0.05(D), 0.10(x), 0.25(+), 0.50(*) have 
been obtained from computer simulations; the lines connect- 
ing these data points merely serve as the guide to the eye.. 




Density 

FIG. 7: Variation of the average speed of the ants in the ant- 
trail model with parallel updating. Same symbols in figs|S| 
and Q correspond to the same values of the parameter /. 

In order to develope a quantitative theory for the flow 
properties of the ant-trail model, we have analyzed the 
spatial organization of the ants by computer simulations. 
Analyzing these observations we concluded that in the 
anomalous regime, loose clusters of ants dominate; the 
term "loose" means that there are small gaps in be- 
tween successive ants in the cluster although the cluster 
appears to be an usual compact cluster if seen from a 
distance. As shown in figlBl the fundamental diagram 
we calculated within the "loose" -cluster approximation 
(LCA) is in good quantitative agreement with the corre- 
sponding data we obtained from computer simulations of 
the model. 



IV. INTRA-CELLULAR TRAFFIC OF 
MOLECULAR MOTORS 

Molecular motors are protein molecules that convert 
the chemical energy, released by the hydrolysis of ATP, 
into mechanical energy required for its forward move- 
ment during intra-cellular transport of vesicular cargo 



FIG. 8: Fundamental diagram (/ = 0.005) of the ant-trail 
model in the LCA (solid curve) s compared with the simula- 
tion data (broken curve). 



|23| . The minimal models developed for explaining the 
mechanism of directed motion of isolated motor proteins 
are based on Brownian ratchets [24|. In such models, 
each motor is represented by a particle. The essential 
features of the detailed mechano-chemistry of the molec- 
ular motor is captured in the Brownian ratchet models 
by a stochastic sequence of successive attachments and 
detachments of the motor with the cytoskeletal filamen- 
tary track (e,g, microtubule in the case of kinesin and 
dynein motors). In the simplest versions of these mod- 
els |25j, in the attached state, the particle representing 
a motor is subjected to a potential that is spatially pe- 
riodic, but each period of which is asymmetric. In the 
detached state the particle executes an unbiased diffusive 
motion. In spite of its simplicity, such a minimal model 
can account for the directed, albeit noisy, movement of 
individual isolated motors. 

To our knowledge, the question of the effects of in- 
teractions of the motors on the intra-cellular traffic was 
addressed theoretically for the first time only a few years 
ago |26| . In that work, the filamentary track was dis- 
cretized in the spirit of the particle-hopping models de- 
scribed above and the motors were represented by field- 
driven particles. Both forward and backward movement 
of the particles were possible and the hopping probability 
of every particle was computed from the local potential. 
Thus, this model was a generalization of ASEP rather 
than TASEP where the hopping probabilities were ob- 
tained from the local potential which itself was time- 
dependent. The fundamental diagram of that model, 
computed imposing periodic boundary conditions, is very 
similar to those shown in fig|2 This observation indi- 
cates that further simplification in the model proposed 
in ref.fHI is possible to develope a minimal model for 
interacting molecular motors. 

Recently, Parmeggiani et al.[8| have, indeed, devel- 
oped such a minimal model for interacting molecular mo- 
tors involved in intra-cellular transport by extending the 
TASEP with open boundary conditions. In this model, 
the molecular motors (e.g., kinesin or dynein) are rep- 
resented by particles whereas the sites for the binding 
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FIG. 9: Schematic representation of the model |g of intra- 
cellular traffic of molecular motors carrying vesicular cargo. 
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FIG. 10: The phase diagram of the model 8| of intra-cellular 
traffic of molecular motors carrying vesicular cargo (repro- 
duced, with permission, from ref.||. The inset shows the 
dependence of the domain wall amplitude on a for different 
values of f3. 
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of the motors with the cytoskeletal tracks (e.g., micro- 
tubules) are represented by a one-dimensional discrete 
lattice. Just as in TASEP, the motors are allowed to hop 
forward, with probability q, provided the site in front is 
empty. However, unlike TASEP, the particles can also 
get "attached" to an empty lattice site, with probability 
A, and "detached" from an occupied site, with probabil- 
ity D (see figini from any site except the end points. The 
state of the system was updated in a random-sequential 
manner. 

To my knowledge, this is the first application of 
TASEP to intra-cellular transport phenomena although 
it is not the first application of TASEP in the domain of 
biological systems; for example, a TASEP-likc model was 
considered earlier for protein synthesis [27l | 

Carrying out Monte-Carlo simulations Parmeggiani et 
al.Q demonstrated a novel phase where low and high 
density regimes, separated from each other by domain 
walls, coexist (see fig llO|) . Using a MFT, they interpreted 
this spatial organization as traffic jam of molecular mo- 
tors. 



FIG. 11: Schematic representation of the different types of 
randomness in particle-hopping models. In (a) the random- 
ness is associated with the track; the hopping probability q 
at the bottleneck (partially hatched region) is smaller than 
the normal hopping probability Q, In (b) the randomness 
is associated with the particles; q\ and §2 being the time- 
independent hopping probabilities of the particles i and j, 
respectively. In (c) the randomness arises from the coupling 
of the dynamics of the hopping particles (filled circle) with 
another non-conserved dynamical variable; the two possible 
states of the non-conserved variable are represented by open 
and filled squares. 



V. DEFECTS AND DISORDER IN 
PARTICLE-HOPPING MODELS 

At least three different types of defects and quenched 
randomness have been considered so far in the context of 
the models of interacting particles driven far from equi- 
librium, (a) First, the randomness may be associated 
with the track on which the particles move; typical ex- 
amples are the bottlenecks on the roads (in the context 
of vehicular traffic) or defects on the microtubules (in 
intra-cellulat transport), etc. For example, as shown in 
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fig lllf a), normal hopping probability at unblocked sites 
is Q whereas that at the bottleneck is q (q < Q). This 
type of quenched defect and disorder of the track leads 
to interesting phase-segregation phenomena (see [| for a 
review) . 

(b) The second type of randomness is associated with 
the hopping particles, rather than with the track. For 
example, the normal hopping probabilities of the par- 
ticles may vary randomly from one particle to another 
(see fig lTTT b)): the hopping probabilities are, however, 
"quenched" random variables, i.e., independent of time. 
In this case, the system is known to be exhibit coarsen- 
ing of queues of the particles and the phenomenon has 
some formal similarities with Bose-Einstcin condensation 
(reviewed in Q). 

Note that in case of the randomness of type (a), the 
hopping probability depends only on the spatial location 
on the track, independent of the identity of the hopping 
particle. On the other hand, in the case of randomness of 
type (b), the hopping probability depends on the hopping 
particle, irrespective of its spatial location on the track. 
In contrast to these two types of randomness, the ran- 
domness in the hopping probabilities of the particles in 
some models arises from the coupling of their dynamics 
with that of another non-conserved dynamical variable. 
For example, in the ant-trail model, the hopping prob- 
ability of an ant depends on the presence or absence of 
pheromone in front of it (see fig lllf c)). Therefore, in 
such models with periodic boundary conditions, a given 
particle may hop from the same site, at different times, 
with different hopping probabilities. 

Defects of either (a) the cytoskeletal filaments or (b) 
the motor proteins or (c) the mechano-chemical coupling 
can cause malfunctiong of the intra-ccllular transport 
leading to various types of diseases [28|. In order to get 
deep insight into the physical origin of such diseases, the 
recent model developed by Parmeggiani et al. [j| has 
been extended 29]. This modeling strategy has opened 
up a new horizon for further unconventional applications 
of non-equilibrium statistical mechanics far beyond the 
traditional borderlines of physics. 

VI. CONCLUSION 

In this paper we began with a brief introduction to 
TASEP [ljj, Ull which is, perhaps, the simplest model 



of systems of interacting driven particles. TASEP, when 
updated in parallel, may be regarded as a special case, 
corresponding to the maximum allowed speed V max = 1 , 
of the NS model [l4| , the minimal model of vehicular traf- 
fic on single- lane highways. We have summarized some of 
the main known results on the fundamental diagram for 
the TASEP with periodic boundary conditions and the 
phase diagram of the TASEP with open boundary condi- 
tions. Then, we have shown how these results for TASEP 
get qualitatively modified by the generalizations or exten- 
sions required to model ant-traffic on ant-trails 0, B| and 
molecular motor traffic on cytoskeletal filaments |8j. 



In the context of the ant-traffic on ant-trails, we have 
established how a combination of analytical and numer- 
ical calculations 0, 0] can account for the unusual shape 
of the fundamental diagram observed in the computer 
simulations of the ant-trail model. We have also pre- 
sented the phase diagram obtained by Parmeggiani et 
al.@ from studies of their recent model for molecular 
motor traffic. This phase diagram suggests the possi- 
bility of coexistence of high-density regions (traffic jam) 
and low-density regions (freely flowing traffic) , separated 
from each other by domain walls, in a novel phase. Fi- 
nally, we have mentioned some ongoing investigations on 
the effects of defects and disorder on molecular motor 
traffic 29] . This trend of research indicates the possibil- 
ity of further unconventional, but very useful, applica- 
tions of statistical physics in biological systems. 
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